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ABSTRACT 


This  thesis  studies  the  effective  thermal  conductivity  of  randomly  oriented, 
percolated  carbon  nanotubes.  To  that  end,  a  multiscale  analysis  approach  was  adopted. 
At  the  nanoscale,  molecular  dynamics  simulation  was  performed  to  determine  the  thermal 
conductivity  coefficient  of  a  single  carbon  nanotube.  Then,  thermal  conductivity  of  two 
carbon  nanotubes  positioned  at  different  angles  were  studied  after  determining  the 
equilibrium  positions  of  the  two  nanotubes  at  various  relative  positions.  Finally,  using 
the  data  obtained  in  the  previous  analyses,  the  effective  thermal  conductivity  of  randomly 
oriented  carbon  nanotubes  was  calculated  using  the  finite  element  model  where  each 
nanotube  was  modeled  as  a  continuous  rod. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Since  being  discovered  by  Iijima  in  1991  [1],  carbon  nanotubes  (CNT’s)  have 
attracted  worldwide  attention  in  both  research  and  industrial  applications.  Carbon 
nanotubes  are  long,  thin  cylinders  of  carbon  macromolecules.  They  can  be  thought  of  as 
a  sheet  of  graphite  (hexagonal  lattice  of  carbon)  rolled  into  a  cylinder  with  a  typical 
diameter  on  the  order  of  0. 14nm.  The  length  of  a  single  walled  nanotube  (SWNT)  can  be 
relatively  large,  up  to  several  microns  long.  The  strong  sp2  bonds  of  carbon  nanotubes 
suggest  an  unusually  high  thermal  conductance,  expected  to  surpass  that  of 
monocrystalline  diamond  [2].  Unfortunately  experimental  attempts  to  quantify  the  heat 
transfer  properties  of  SWNT  have  encountered  many  difficulties.  However  molecular 
dynamics  (MD)  simulations  has  proven  to  be  a  useful  contributor  to  the  ongoing  catalog 
of  CNT  heat  transfer  properties. 

The  thennal  conductivity  is  a  thennophysical  property  of  a  material  that  indicates 
its  ability  to  conduct  heat.  The  magnitude  of  thennal  conductivity  varies  over  wide 
ranges  for  different  materials.  Thermal  conductivity  of  a  material  depends  on  its 
chemical  composition,  physical  structure,  and  state.  In  addition,  it  varies  with  the 
pressure  and  temperature  at  which  the  material  is  exposed.  However,  pressure  is  far  less 
influential,  so  the  dependence  is  tabulated  as  a  function  of  temperature.  In  anisotropic 
materials,  such  as  carbon  nanotubes,  the  thennal  conductivity  also  varies  with  direction. 
Given  its  small  size  and  potentially  high  thermal  conductivity,  researchers  have  focused 
on  this  material  as  an  efficient  heat  conductor  in  nanodevices  and  electronics. 

Both  diamond  and  in-plane  graphite  display  very  high  thermal  conductivities  at 
moderate  temperatures;  for  example  4900  and  3870  (W  m-1  K-l)  at  173.2  K, 
respectively  [3],  Both  are  made  of  a  covalently  bonded  network  of  carbon  atoms. 
Having  a  similar  lattice  structure,  it  is  expected  that  carbon  nanotubes  will  also  have  a 
high  thennal  conductivity.  In  fact  many  theoretical  and  experimental  results  suggest  that 
the  thermal  conductivity  of  CNTs  may  be  larger  than  those  of  both  diamond  and  in-plane 
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graphite.  Thus,  an  accurate  evaluation  of  the  thennal  conductivity  of  nanotubes  is  an 
important  step  in  the  development  of  useful  nanoscale  devices. 

Theoretical  calculations  have  predicted  the  thermal  conductivity  of  a  single 
nanotube  to  be  as  high  as  6600  W  m-1  K-l  at  room  temperature  [2],  However 
experimental  measurements  have  not  been  as  high.  Experimental  measurements  of  “mat” 
samples  revealed  a  room  temperature  conductivity  of  approximately  35  W  m-1  K-l  [4]. 
This  degradation  is  due  to  the  disordered  and  highly  tangled  nature  of  the  sample.  In 
addition,  CNT  has  much  lower  thermal  conductivity  in  the  transverse  direction  than  the 
longitudinal  direction.  These  extreme  differences  in  thermal  conductivity  are  indicative 
of  the  fact  that  multiwalled  nanotubes  (MWNT),  bundles  and  ropes  exhibit  varying 
degrees  of  thermal  conductivity. 

Recent  publications  have  reported  both  thennal  and  electrical  properties 
enhancements  of  single  wall  carbon  nanotube-polymer  composites.  Scientific  results 
show  that  these  enhancements  are  a  function  of  temperature,  nanotube  distribution  and 
concentration.  One  report  presented  by  Biercuk  et  al.  [5]  showed  epoxy  samples  loaded 
with  1  wt%  unpurified  SWNT  material  increased  thermal  conductivity  by  70%  at  40K  in 
comparison  to  pristine  epoxy.  The  conductivity  further  increased  to  125%  at  room 
temperature.  Figure  1  illustrates  the  enhancement  in  thermal  conductivity  relative  to 
pristine  epoxy  as  a  function  of  SWNT  and  vapor  grown  carbon  fibers  (VGCF). 
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Figure  1.  Enhancement  in  thermal  conductivity  relative  to  pristine  epoxy  as  a  function 

of  SWNT  and  VGCF  loading  (From  Ref.  [5]) 


Meanwhile,  researchers  at  Florida  State  University  studied  the  effect  of  magnetic 
alignment  on  composite  materials  with  3  wt%  SWNT  loading.  Even  without  magnetic 
field  processing,  the  thermal  conductivity  increased  by  300%.  Magnetic  aligmnent 
further  increased  the  thennal  conductivity  by  additional  10%  [6]. 
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Figure  2.  Thermal  conductivity  of  CNT-epoxy  composites  magnetically  processed  at  0 
and  25  T,  compared  with  the  neat  epoxy  control  sample  (also  processed  at  25  T) 

(From  Ref.  [6]) 

The  above  experimental  observations  provided  significant  evidence  of  the  use  of 
CNTs  as  a  conductor  of  heat  in  composite  materials,  as  well  as  in  nano-scale  devices. 
Thus,  it  is  necessary  to  determine  how  orientation  and  dispersion  affect  heat  conduction 
in  nano  tube  systems.  Fortunately,  computational  techniques  for  modeling  and  simulation 
which  bridge  continuum  mechanics,  quantum  mechanics  and  molecular  dynamics  are 
available  tools  that  predict  and  explain  thermal  properties  of  various  CNTs. 

B.  OBJECTIVES 

In  this  thesis,  we  employ  a  combination  of  the  classical  MD  simulation  method 
and  the  finite  element  method  (FEM)  to  analyze  the  thennal  conductivity  in  the 
longitudinal  and  transverse  directions  of  randomly  dispersed  SWNTs. 

To  this  aim,  a  multiscale  analysis  approach  was  adopted.  At  the  nanoscale, 
molecular  dynamics  simulation  was  performed  to  determine  the  heat  conduction 
coefficient  of  a  single  carbon  nanotube  using  the  Tersoff-Brenner  potential  [7].  Then, 
thennal  conductivity  of  two  carbon  nanotubes  positioned  at  different  angles  were  studied 
after  detennining  the  equilibrium  positions  of  the  two  nanotubes  at  various  relative 
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positions  using  the  Lennard- Jones  pair  potential  for  van  der  Waals  interaction  between 
adjacent  SWNTs.  Finally,  using  the  data  obtained  in  the  previous  analyses,  the  effective 
thennal  conductivity  of  randomly  oriented  carbon  nanotubes  was  calculated  using  the 
finite  element  method  where  each  nanotube  was  modeled  as  a  continuous  rod.  As  a  final 
step,  an  estimation  of  the  validity  of  the  molecular  dynamics  and  FEM  applied  in  this 
thesis  is  made  by  comparative  analysis  of  the  results  with  experimental  data. 
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II.  MODELING,  RELATED  THEORY,  AND  METHODOLOGY 


The  properties  exhibited  by  carbon  nanotubes  are  highly  dependent  upon  the  basic 
geometrical  structure  of  the  molecule.  Understanding  of  this  dependence  is  an  essential 
step  in  application  of  the  MD  simulation  of  a  given  CNT  system.  For  example,  the  CNT 
structure  should  be  analyzed  in  its  lowest  energy  state.  This  chapter  explains  the 
synthesis  of  SWNTs,  describes  the  basic  structural  properties  of  carbon  nanotubes,  and 
then  presents  the  resultant  initial  geometrical  models  of  SWNT  which  allow  for  the 
application  of  molecular  dynamics  and  finite  element  modeling  approach. 

A.  LABORATORY  SYNTHESIS  OF  CARBON  NANOTUBES 

There  are  three  common  techniques  used  to  produce  carbon  nanotubes  of  both  the 
single-wall  and  multi-wall  varieties:  (1)  arc  discharge  (2)  laser  ablation  and  (3)  chemical 
vapor  deposition.  Purification  of  the  tubes  can  be  achieved  by  a  variety  of  techniques: 
oxidation,  acid  treatment,  annealing,  sonication,  and  filtering  and  functionlisation 
techniques.  There  are  pros  and  cons  to  be  considered  for  each  of  these  techniques. 
However  the  ultimate  goal  is  to  create  an  economically  feasible  large  scale  production 
technique  that  yields  replicable  pure  results.  Pure  nano  tubes  are  essential  to  the 
development  and  functionality  of  nano-devices  and  composite  materials.  It  is  widely 
accepted  in  the  nanotechnology  field  that  impure,  flawed  nanotubes  have  significantly 
degraded  physical  properties  in  comparison  to  pure  nanotubes. 

The  arc  discharge  method  is  the  most  common  and  simplest  way  to  produce 
carbon  nanotubes.  However  the  product  consists  of  a  mixture  of  components  and 
requires  separating  nanotubes  from  the  byproducts  of  the  reaction.  This  method  creates 
nanotubes  through  arc-vaporization  of  two  carbon  rods  placed  end  to  end,  separated  by 
approximately  1mm,  in  an  enclosure  that  is  usually  filled  with  inert  gas  (helium,  argon)  at 
low  pressure  (between  50  and  700  mbar).  A  direct  current  of  50  to  100  A  driven  by 
approximately  20  V  creates  a  high  temperature  discharge  between  the  two  electrodes. 
The  discharge  vaporizes  one  of  the  carbon  rods  and  forms  a  small  rod  shaped  deposit  on 
the  other  rod.  Producing  nanotubes  in  high  yield  depends  on  the  uniformity  of  the  plasma 
arc  and  the  temperature  of  the  deposit  form  on  the  carbon  electrode. 
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The  laser  ablation  method  is  similar  to  the  arc  discharge  method  in  that  the 
reaction  takes  place  in  the  presence  of  an  inert  gas  and  a  catalyst.  In  this  method, 
graphite  electrodes  are  vaporized  by  a  laser  in  an  oven  of  1200°C.  The  oven  is  filled  with 
the  inert  gas  to  limit  the  pressure  of  the  reaction.  As  the  laser  makes  contact  with 
graphite,  a  hot  vapor  plume  forms,  expands  and  cools  rapidly.  As  the  vaporized  particles 
cool,  small  carbon  molecules  and  atoms  condense  to  fonn  a  mixture  of  carbon  nanotubes 
and  carbon  nanoparticles.  The  use  of  pure  graphite  electrodes  yields  MWNTs;  and 
SWNTs  can  be  synthesized  from  a  mixture  of  graphite  and  Co,  Ni,  or  Fe.  Laser 
vaporization  results  in  a  higher  purer  yield  of  SWNTs  than  those  produced  by  the  arc 
discharge  method. 

The  third  method,  chemical  vapor  deposition  (CVD)  synthesis  requires  a  gaseous 
form  of  carbon  such  as  methane  or  carbon  monoxide  and  an  energy  source,  such  as 
plasma.  The  energy  source  is  used  to  separate  the  carbon  molecule  into  reactive  atomic 
carbon.  The  carbon  then  diffuses  towards  a  pre-treated  substrate,  where  the  atoms 
combine  to  form  the  carbon  nano  tubes.  Growth  of  SWNTs  or  MWNTs  is  controlled  by 
applying  an  appropriate  catalyst  to  the  substrate. 

The  exact  mechanism  by  which  CNT’s  are  formed  is  not  known.  However  the 
three  synthesis  techniques  all  require  a  carbon  based  precursor,  a  catalysis,  and  intense, 
high  energy  reaction  conditions. 

B.  MD  THEORY  AND  TERSOFF-BRENNER  (T-B)  POTENTIAL 

1.  Validity  of  the  Classical  MD  Simulation  Method 

In  molecular  dynamics  (MD)  simulations,  atoms  are  treated  as  discrete  particles 
whose  trajectories  are  calculated  by  numerically  integrating  classical  equations  of  motion 
subject  to  some  given  interatomic  forces.  The  interatomic  interactions  are  calculated 
from  one  of  two  methods.  The  first  approach  requires  a  summation  of  nuclear  repulsions 
combined  with  electronic  interactions  determined  from  first  principles.  This  method  is 
commonly  referred  to  as  a  quantum  mechanical  MD  because  it  considers  the  influence 
electrons  have  on  physical  properties.  Quantum  MD  simulations  are  computationally 
more  difficult  and  time  consuming  because  of  the  inclusion  of  electron  interactions.  In 
the  second  approach,  quantum  mechanical  electrons  are  ignored.  The  energy  and 

interatomic  forces  are  a  function  of  atomic  position  only.  This  approach  is  known  as 
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classical  MD,  and  is  suitable  for  analyzing  the  mechanical  properties  of  molecules  when 
quantum  effects  are  not  significant.  Under  the  appropriate  assumptions  energetically 
stable  structures  like  carbon  nanotubes  can  be  modeled  with  the  classical  approach  using 
Newton’s  second  law,  the  Hamiltonian  equation  of  motion,  and  an  appropriate  potential 
energy  function.  In  this  thesis,  the  classical  MD  simulation  method  was  employed  to 
calculate  the  thermal  conductivity  of  a  SWNT. 

2.  Newton-Hamiltonian  Dynamics  for  the  Classical  MD  Simulation 

In  the  classical  molecular  dynamics  method,  the  equations  of  motion  (Newton’s 
2nd  Law)  are  solved  for  atoms  as: 


F.  =  mi]  =  m 


(1) 


Here  r,  and  Ft  are  position  vector,  and  force  vector  of  atom  i.  The  position  vector  locates 
the  atom  with  respect  to  the  origin  of  the  coordinate  system.  The  mass  is  assumed  to  be 
independent  of  position  velocity,  and  time. 

For  an  isolated  system,  the  total  energy  E  is  conserved.  Therefore,  we  can 
identify  the  total  energy  as  the  Hamiltonian  H.  For  N-atoms,  H  takes  the  form: 


1  N 

H(ri  > Pi  )=^X^2  +u(ri  )=E 


(2) 


where  pt  is  the  momentum  of  an  atom  and  the  potential  energy  U  results  from  interatomic 
interactions.  Taking  the  derivative  of  equation  (2)  with  respect  to  the  position  r,  the 
explicit  relationship  between  the  Hamiltonian  and  the  potential  energy  is: 

dH  _dU 

dr  dr  (3) 

Therefore,  the  resultant  Hamilton’s  equations  of  motion  are: 

dll  _  pt  _  . 

dpt  m  '  (4) 
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dH 

dr 


(5) 


=  -Pi  =  ~mn 


Substituting  equation  (3)  into  equation  (5)  Newton’s  law  becomes: 

dH  _  8U 
dr  dr 


(6) 


Equations  (6)  shows  any  conservative  force  can  be  written  as  the  negative  gradient  of 
some  potential  energy  function  U,  and  the  forces  from  all  of  the  other  atoms  determine 
the  accelerations  of  each  atom  in  the  given  system  [8]. 

3.  Tersoff-Brenner  (T-B)  Type  Potential 

To  model  the  bonded  interaction  among  the  carbon  atoms  of  the  nanotube,  the 
classical  Tersoff-Brenner  potential  was  used  [7].  The  Tersoff-Brenner  potential  is 
specially  suited  for  carbon-based  systems  and  has  been  used  in  a  wide  variety  of 
applications  yielding  results  in  agreement  with  experimental  observations  [9], 

A  universal  analytic  many-body  force  field  function  that  works  for  all  materials  in 
all  scenarios  does  not  currently  exist.  However,  the  Tersoff-Brenner  potential  for  carbon- 
based  systems  allows  for  reactive  short-range  bonded  interactions.  Specifically,  chemical 
bonds  can  form  and  break  during  the  course  of  the  simulation.  This  significant  advantage 
was  contributed  by  Brenner  when  he  introduced  a  bond  function  to  Tersoff  s  original 
analytic  potential  energy  function  [10].  The  resulting  empirical  equation  is  a  pair  wise 
exponential  potential  function  of  the  atomic  separation,  angle,  and  the  number  of  bonds 
given  by: 


^  i 


where  the  repulsive  and  attractive  pair  terms.  Vr  and  VA  are  given  by: 

VK{')j)  =  fArij) DT  /(<Sv  - ’ exp (-^25^ (r - Rf ) 

VA  (rv)  =  fj  (rv)D<,%  /OS,  - 1) •  exp ( -^2AST p.. (r  - Rf  ) 


(7) 

(8) 

(9) 

(10) 
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D,  S,  /?,  R  are  empirical  parameters  given  in  Table  1,  and  the  function  R  (r) ,  restricts  the 
pair  potential  to  nearest  neighboring  atoms,  given  by: 


Mr)  =  \ 


1, 


f  /r(r-R(1})^ 

1  +  cos 

V  v  J 

r(2)-rw 

l  v  11  ) 

/  2 


0, 


r  <  R, 


(i) 


>  SP<r<$ 


(2) 


r  <R 


(2) 


(ID 


Bij  represents  a  many-body  coupling  between  the  bond  from  atom  i  to  atom  j  and  the 
local  environment  of  atom  i. 


= 


~\-Si 


»] 


(12) 


G(9ijk)  =  ci(i 


i+V 

dl 


d\  +(l  +  cos^)‘ 


(13) 


G  is  a  function  of  6,  the  angle  between  atoms  i  -  j  and  i  -  k  bonds  and  other  constant 
variables.  These  equations  require  a  total  of  1 1  fitting  parameters  for  carbon  atoms.  The 
values  assigned  to  each  of  these  parameters  are  given  in  Table  1,  where  the  subscripts 
“cc”  means  that  those  parameters  are  used  for  just  carbon-carbon  bonding  network 
without  other  elements  such  as  hydrogen. 
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Table  1.  Optimized  Parameters  for  the  T-B  Type  Potential.  [7] 


Parameter 

Value 

Parameter 

Value 

Parameter 

Value 

1.42  A 

<*«c 

0.0 

dl 

3.52 

D? 

6.0  eV 

1.7  A 

Pcc 

2.1  A 

*2? 

2.0  A 

1.22 

ao 

0.00020813 

&cc 

0.5 

co 

3302 

The  T-B  potential  energy  function  of  two  carbon  atoms  can  be  plotted  as  a 
continuous  function  from  equations  (9)  through  (15).  The  effective  carbon  bond  length, 
r0~1.42A  is  the  equilibrium  separation  distance  between  two  carbon  atoms.  The 
equilibrium  of  two  carbon  atoms  corresponds  to  the  minimum  potential  energy  of  the 
system.  In  order  for  any  CNT  structure  to  form,  the  potential  energy  per  each  carbon 
atom  must  attain  a  minimum  energy  value.  This  minimum  energy  is  defined  as  the  bond 
energy  and  is  known  to  be  e~2.5eV  for  the  carbon-carbon  bond  of  CNTs.  This  value  is 
consistent  with  experimental  values  of  C-C  tight  bonding  overlap  energy  [11]. 

The  concept  of  bond  formation  occurring  at  the  lowest  potential  energy  is 
consistent  throughout  the  world  of  chemistry.  A  bounded  system,  such  as  a  carbon 
nanotube,  has  a  lower  potential  energy  than  its  constituent  parts,  which  is  what  keeps  the 
system  together.  The  carbon  atoms  exert  a  force  on  each  other  in  order  to  position 
themselves  in  the  lowest  potential  energy  state,  graphically  shown  as  the  bottom  of  the 
potential  well.  This  force  is  equivalent  to  the  negative  gradient  of  the  potential  energy 
function  U  given  in  equation  (8). 
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4.  Gear’s  Predictor-Corrector  Method  for  Solving  the  EOM 

All  MD  calculations  require  accurate  and  stable  integrators  in  order  to  solve  the 
equations  of  motion.  Integrators  are  evaluated  by  stability  and  the  ability  to  reproduce 
certain  time  and  space  correlations  to  a  sufficient  degree  of  accuracy.  The  Beeman 
algorithm  and  Verlet  leap-frog  method  are  classical  finite  difference  approaches.  The 
velocity  of  each  molecule  is  calculated  first,  followed  by  the  position  of  the  molecule. 
However,  in  molecular  dynamics,  the  motion  of  a  particle  is  over  a  very  large  number  of 
time  steps,  and  the  elaborate  integration  schemes  of  the  Gear’s  Predictor-Corrector 
method  has  proven  to  be  more  accurate  in  each  time  step. 

Gear’s  Predictor-Corrector  algorithms  consist  of  three  steps;  prediction, 
evaluation,  and  correction: 

a.  Prediction 

From  the  positions  and  their  time  derivatives  up  to  a  certain  order  (fifth- 
order  in  this  case),  all  known  at  time  t,  one  "predicts"  the  same  quantities  at  time  t  +At 
by  means  of  a  Taylor  expansion.  Among  these  quantities  are  acceleration,  velocity  and 
position  of  the  carbon  atom. 


ri(t  +  At)  =  ri(t)  +  ri(t)At  +  ri(t)  f 


ri(t  +  At)  =  ri(t)  +  ri(t)At  +  n(t)(  J  +r(  )(tj 

(A,)4 

4! 

(15) 

ri(t  +  At)  =  ri(t)  +  ri(t)At  +  ri(  )  (t)^  J 

(16) 

' r(t  +  At)  =  r{t)  +  rf4)  (t)  At  +  rt f 1  ( 

,)(aT 
'  ’  2! 

(17) 

r^\t  +  At)  =  rf-A)  (t)  +  r;.(5)  (t)At 

(18) 

#;(5)  (f  + AO  =  r/5)(f)  (19) 


13 


b.  Evaluation 

The  interatomic  force  F,  on  each  atom  is  computed  taking  the  gradient  of 
the  potential  at  the  predicted  positions.  The  force  on  each  atom  is  given  by  equation  (20), 
where  U(ry)  is  the  continuous  potential  energy  function  that  acts  between  atoms  i  and  j. 
The  resulting  acceleration  will  be  in  general  different  from  the  "predicted  acceleration". 
The  difference  between  the  two  constitutes  an  "error  signal”. 


ySUtjjX 


(20) 


c.  Correction 

This  error  signal  A r.  is  used  to  "correct"  positions  and  their  derivatives. 


All  the  corrections  are  proportional  to  the  error  signal,  the  coefficient  of  proportionality 
being  a  "magic  number"  determined  to  maximize  the  stability  of  the  algorithm  [12]. 


A r.  =  r.  {t  +  At)  -  r.p  [t  +  At) 


(21) 


5.  Simulation  Time  Step 

By  iterating  the  three  steps  of  Gear’s  method  at  every  time  step,  the  MD 
simulation  is  propagated  through  time  at  intervals  of  At.  Using  the  relationship  between 
bond  length  r0  and  the  bond  energy  s,  a  reasonable  simulation  time  step  At  of  0.01  ps  was 
estimated  for  the  MD  simulation.  The  methods  required  for  this  estimation  are  outlined 
below. 

The  kinetic  energy  of  an  atom  in  an  isolated  system  should  be  approximately 
equal  to  the  potential  energy  of  the  atom.  Assuming  that  one  atom  can  move  within  the 
bond  length  ra  during  the  time  step,  then  the  kinetic  energy  and  the  velocity  v  of  the  atom 
is  written  as: 


E 


k 


(22) 


At 

From  equations  (22)  and  (23),  the  maximum  possible  time  step  value  is  given  by: 

At  =  r0y/m/2e  <  0.02  ps 
14 


(23) 


(24) 


where  m  is  the  mass  of  the  carbon  atom,  and  rQ  and  a  the  carbon  bond  length  and  C-C 
tight  bonding  overlap  energy.  Equation  (24)  confirms  that  At=  O.Olps  is  a  reasonable 
time  step  for  the  MD  code. 

C.  EQUILIBRIUM  MOLECULAR  DYNAMICS  USING  AUTO¬ 
CORRELATION  LORMULA  (GREEN-KUBO  FORMULA) 

This  section  describes  the  method  of  calculating  thermal  conductivity  of  a  SWNT 

by  introducing  a  heat  flux  autocorrelation  function.  Fourier’s  law,  which  describes 

macroscopic  thermal  conduction,  is  not  appropriate  for  low-dimensional  nano  systems. 

Fourier’s  law  states  that  the  heat  flux  j  is  related  to  the  temperature  gradient  as  j  =  -kV  T, 

where  k  is  the  thennal  conductivity  tensor  and  T  is  the  temperature  distribution. 

Therefore  a  simple  approach  to  studying  thennal  conduction  in  SWNTs  is  to  apply  a 

temperature  gradient,  measure  heat  flux  along  the  axial  direction,  then  calculate  the 

thennal  conductance.  However  the  linearity  between  heat  flux  and  temperature  gradient 

as  defined  by  Fourier’s  law  does  not  hold  due  to  the  small  dimensions  of  a  CNT. 

Therefore  the  Green-Kubo  fluctuation-dissipation  theorem  [13]  is  employed  in  this  work. 

In  an  equilibrium  system,  Green-Kubo  formula  takes  the  form: 

(25) 

where  k#  is  the  Boltzmann  constant,  V  is  the  volume,  T  is  the  temperature  of  the  sample 
and  the  angular  brackets  denote  an  ensemble  average.  The  microscopic  heat  flux  vector 
J(t)  is  defined  by: 

=  (26) 

at  i 

where  r/t)  is  the  time-dependent  coordinate  of  atom  i.  According  to  the  MD  theory,  the 
total  potential  energy  can  be  expressed  as  a  sum  of  binding  energies  of  individual  atoms, 
therefore  the  site  energy  e/t)  can  be  taken  to  be: 

A  =^A,2  +  ^ZM(f/)  (2?) 

In  the  above  equation,  v,  represents  the  velocity  vector  and  u(ry)  is  the  Tersoff- 

Brenner  potential  function  described  previously. 
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Generally  two  physical  mechanisms  contribute  to  thermal  conduction  in  CNTs:  (i) 
electron-phonon  interactions,  which  depend  on  electronic  band  structures  and  the  electron 
scattering  process  and  (ii)  phonon-phonon  interactions,  which  depend  on  the  vibrational 
modes  of  the  lattice  [14].  For  CNTs,  the  electronic  contribution  to  the  thermal 
conductivity  is  negligible  due  to  the  low  density  of  free  charge  carriers  [15].  Phonon- 
phonon  interactions  dominate  the  overall  thermal  conductivity  and  can  be  studied  by 
classical  MD 

D.  LENNARD-JONES  POTENTIAL 

Non-bonded  intermolecular  interaction  of  carbon  atoms  in  SWNT’s  can  be 
represented  by  the  Lennard-Jones  potential.  The  Lennard-Jones  potential  is  a  simple 
mathematical  model  that  represents  two  distinct  forces  of  neutral  atoms  or  molecules. 
When  the  separation  r  between  the  two  molecules  is  small,  Pauli  repulsive  forces 
dominate  resulting  in  strongly  positive  potential.  In  contrast,  when  the  separation  r 
increases  an  attractive  force  known  as  van  der  Waals  dominates.  The  general  equation  is 
expressed  as: 


V(r)  =  4s 


f 

12 

f  ( 7 

6 

yr  j 

\r  J 

where,  r  is  the  distance  between  interacting  atoms,  s  and  a  are  Lennard-Jones  parameters 
for  energy  and  equilibrium  distance  respectively.  For  carbon  atoms  the  Lennard-Jones 
parameters  are  s  =  0.00778  kcal/mol  and  a  =3.4  A  [16],  The  Lennard-Jones  model 
provides  a  smooth  transition  between  (a)  the  attractive  van  der  Waals  force  of  an 
approaching  pair  of  atoms  from  a  certain  distance  and  (b)  the  repulsive  force  when  the 
distance  of  the  interacting  atoms  becomes  less  that  the  sum  of  their  contact  radii.  At 
equilibrium,  the  pair  of  atoms  or  molecules  tends  to  go  toward  a  separation  corresponding 
to  the  minimum  of  the  Lennard-Jones  potential. 
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III.  SIMULATION 


A.  MODELING  OF  SWNT  AND  BASIC  STRUCTURES 

1.  Basic  Structures  of  Carbon  Nanotubes 

Carbon  nanotubes  are  long,  thin  cylindrical  shell  structures  based  on  the 
hexagonal  lattice  structure  of  carbon  allotrope,  graphite.  The  length  to  diameter  ratio  can 
be  extremely  large,  and  the  ends  are  “capped”  with  half-dome  shaped  fullerene 
molecules. 

The  single-walled  carbon  nanotube  is  considered  to  be  the  cornerstone  of  the 
nano  tube  family.  Other  known  variants  are  multi- walled  nano  tubes  (MWNT)  and 
ordered  SWNT  arrays  known  as  nanoropes.  Most  of  the  existing  experimental  methods 
have  been  successful  at  yielding  MWNTs,  however  isolation  of  pure  SWNTs  have  been 
the  focus  of  many  recent  experimental  techniques.  If  an  adequate  understanding  of 
SWNTs  is  obtained,  the  cataloging  of  other  variant  properties  will  be  much  easier.  The 
following  outline  will  concentrate  on  the  basic  structural  characteristics  of  SWNTs. 


Figure  3.  Graphite  layer  with  carbon  atoms  labeled  using  (n,m)  notation[6] 

To  understand  the  structural  features  of  SWNT,  Figure  3  shows  a  two 
dimensional  honeycomb  structure  of  a  graphene  sheet.  The  vectors  shown  represent  the 
chiral  vectors  of  a  nanotube  (C).  A  chiral  vector  is  defined  as,  C  =  nd]  +ma2 ,  where  a{ 
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and  a-,  are  unit  vectors.  The  vector  represents  two  atoms  in  the  graphene  sheet,  one  of 

which  serves  as  the  origin.  If  the  sheet  is  rolled  until  the  two  atoms  coincide,  the  vector 
pointing  from  the  first  atom  towards  the  other  is  the  chiral  vector  and  its  length  is  equal  to 
the  circumference  of  the  SWNT.  Depending  upon  the  combination  of  the  indices  (n,  in), 
three  types  of  CNTs  is  possible.  An  (m,  m)  combination  forms  the  “armchair”  tube,  (n,0) 
is  known  as  “zig-zag”,  and  (n,  m)  refers  to  a  chiral  nanotube.  The  pseudonym  illustrates 
the  shape  of  the  nanotube  perpendicular  to  the  tube  axis 

The  physical  properties  of  carbon  nanotubes  are  largely  determined  by  the  chiral 
indices  (n,  in).  The  diameter  of  nanotubes,  ^expressed  as  a  function  of  n  and  m  is  given 
by  the  following  equation. 


d,  = 


0.246 


n 


+  n  ■  m  +  m1)  ( nm ) 


(28) 


The  chiral  angle  is  the  angle  between  C  and  ax  ,  given  by  equation  (29). 


6  =  sin  1  \ 


31/2  • 


m 


2( 


2  2 \V2 

n  +  n  ■  m  +  in 


(29) 


In  addition,  a  SWNT  can  be  either  a  metal,  semiconductor  or  small-gap  semiconductor 
depending  on  the  (n,  in)  combination.  For  instance,  nanotubes  with  m  -  n  ^  3  x  integer 
results  in  a  semi-conducting  nanotube.  Nanotubes  with  m  -  n  =  3  x  integer,  is  found  to  be 
metallic[17].  Although  these  characterizes  are  very  important  to  the  study  of  electronic 
properties  of  SWNT,  only  a  weak  dependence  of  thennal  conductivity  can  be  contributed 
to  chirality [  1 8] . 

2.  SWNT  Model 

Three  models  of  SWNTs  were  created  in  order  to  analyze  the  heat  transfer 
properties.  A  single  (6,  6)  armchair  CNT  was  modeled  in  order  to  determine  the  thermal 
conductivity  of  a  typical  SWNT.  A  second  CNT,  of  equivalent  dimensions  was  then 
created  to  analysis  the  effects  of  intennolecular  forces  on  thermal  conductivity.  Figure  4 
shows  the  side  view  of  two  parallel  SWNTs  each  consisting  of  366  carbon  atoms,  which 
are  represented  by  individual  circles  through  the  visualization  MATLAB  algorithm. 
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Figure  4.  Geometric  model  two  parallel  (6,  6)  armchair  SWNT’s 

The  C-C  distance  used  in  this  model  is  1.421  A,  and  the  radius  of  the  tube  is 
4.07A.  The  length  to  diameter  ratio  was  limited  to  5: 1  in  order  to  decrease  computational 
costs.  A  third  model  was  created  to  determine  the  effective  thermal  conductivity  of 
randomly  dispersed  SWNTs.  The  CNT’s  are  represented  as  1-D  rods  randomly  oriented 
in  unit  square.  Figure  5  shows  a  typical  dispersion  of  50  SWNT  with  nondimensional 
length  of  0.5. 
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Figure  5.  1-D  representation  of  50  randomly  dispersed  SWNTs 


B.  SIMULATION  PROCEDURES 

All  programming  language  required  for  the  computer  simulations  in  this  study 
were  written  using  Mathworks-MATLAB.  MATLAB  is  a  high-level  technical 
computing  language  and  interactive  environment  for  algorithm  development,  data 
visualization,  data  analysis,  and  numeric  computation. 

1.  Thermal  Conductivity  Coefficient  of  a  Single  Carbon  Nanotube 

An  armchair  (6,  6)  SWNT  structure  was  modeled,  having  a  diameter  of  8.14  A 
and  length  to  diameter  ratio  of  5: 1 .  The  model  consisted  of  366  atoms. 
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Armchair  (6,6)  SWNT,  No.  Atoms  =  294 


Figure  6.  Geometric  model  (6,  6)  armchair  SWNT 

In  this  molecular  dynamics  simulation  program,  the  Tersoff  Brenner  potential  was 
employed  as  the  potential  function  with  the  parameters  set  to  those  outlined  in  Table  1. 
The  Gear’s  Predictor-Corrector  method  was  adapted  to  integrate  the  equation  of  motion 
with  each  run  consisting  of  10000  time  steps  of  0.01  ps.  Periodic  boundary  conditions 
were  used  to  model  the  microscopic  system.  The  thermal  conductivity  coefficient  in  both 
the  axial  and  radial  direction  was  found  from  the  Green-Kubo  autocorrelation  function 
outlined  in  equation  (25). 

2.  Thermal  Conductivity  of  Two  Carbon  Nanotubes 

Prior  to  determining  the  thermal  conductivity  of  two  adjacent  SWNTs,  the  non- 
bonded  equilibrium  distance  was  determined  from  the  Lennard-Jones  potential  function. 
For  this  simulation  two  equivalent  armchair  (6,  6)  SWNTs  were  modeled.  The 
equilibrium  position  was  detennined  by  incrementally  increasing  the  relative  distance 
between  the  two  CNTs  from  0  to  10  while  simultaneously  calculating  the  LJ  potential  for 
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the  system.  These  measurements  were  calculated  at  fixed  relative  angles  of  15,  25,  30, 
45,  60,  75  and  90  degrees.  The  equilibrium  position  was  then  used  in  the  detennination 
of  the  thermal  conductivity  of  two  adjacent  SWNTs.  The  molecular  dynamic  methods 
described  in  section  one  where  simply  repeated  for  two  adjacent  CNTs.  By  doing  so,  the 
affect  of  non-bonded  carbon-carbon  interactions  is  represented  in  the  longitudinal  and 
transverse  thermal  conductivities  calculations. 

3.  Effective  Thermal  Conductivity 

The  effective  thennal  conductivity  of  randomly  oriented  carbon  nanotubes  was 
calculated  using  finite  element  analysis.  Unlike  macroscopic  three-dimensional  material, 
carbon  nanotubes  can  be  thought  of  as  quasi  one-dimensional  rods.  Therefore,  to  simulate 
thennal  contact  between  adjacent  SWNTs  each  nanotube  was  modeled  as  a  continuous 
rod.  The  total  number  of  tubes  was  varied  between  25  and  100  in  increments  of  25. 
Additionally,  the  non-dimensional  length  of  the  CNT  was  set  at  0.5,  0.6,  0.7,  and  0.8  to 
determine  the  effect  length  has  on  thermal  conductivity.  Finally,  the  non-dimensional 
lengths  were  decreased  by  50%  while  doubling  the  number  of  CNTs.  This  procedure 
allowed  for  the  analysis  of  tube  length  while  keeping  the  total  volume  of  CNTs  constant. 
The  thermal  conductivity  values  determined  in  the  first  MD  simulation  is  used  to  describe 
the  material  properties  of  randomly  dispersed  nanotubes  in  this  final  simulation. 

MATLAB’s  “rand”  function  was  used  to  generate  random  position  coordinates  of 
the  initial  tube.  Subsequent  tubes  were  then  added  in  a  random  manner.  Concurrently, 
the  tubes  were  normalized  and  repositioned  to  ensure  that  they  remained  with  the  domain 
of  a  unit  square.  To  ensure  percolation  between  the  temperature  boundary  conditions  and 
to  prevent  the  generation  of  an  isolated  tube,  each  tube  was  required  to  overlap  with 
another  tube  within  the  confines  of  the  unit  square.  Conventional  finite  element 
procedures  were  then  employed  for  the  heat  conduction  analysis. 

For  finite  element  analysis  of  microscale  heat  transfer  of  CNTs,  nodal  coordinates 
were  extracted  from  the  end  point  coordinates  of  the  1-D  rods  as  well  as  from  the 
intersecting  coordinates  of  the  randomly  overlapped  tubes.  Hence,  a  truss  element  with 
two  nodes  was  used  to  represent  the  basic  CNT  structure.  From  Fourier’s  law,  the 
governing  equation  for  one-dimensional  heat  flux  can  be  written: 
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J  =  -kAAT  (30) 


where  k  is  the  thermal  conductivity  tensor  ,A  is  the  cross  section  area,  and  T  is  the 
temperature  distribution.  The  finite  element  formulation  leads  to  a  generic  matrix 
equation  in  the  form  of 

[K]{u}  =  {F}  (31) 

K  being  the  stiffness  matrix,  and  u  the  temperature  vector.  The  heat  flux  vector,  F, 
represents  the  contributions  of  the  heat  source  and  the  boundary  conditions.  Solving  for 
the  nodal  temperatures,  we  are  able  to  calculate  horizontal  and  vertical  heat  flux  of  the 
unit  square.  The  effective  thermal  conductivity  was  found  by  substituting  the  heat  flux 
solutions  into  equation  (30)  and  solving  for  the  thennal  conductivity  coefficient. 
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IV.  RESULTS  AND  DISCUSSION 


A.  THERMAL  TRANSFER  COEFFICIENT  OF  A  SINGLE  CARBON 
NANOTUBE 

In  this  work,  several  annchair  CNTs  are  simulated. 


1 

fable  2.  Thermal  Conductivity  of  SWNTs 

Armchair  Indices 

Longitudinal  Direction 

Reference  Values 

(5,5) 

3325.8  W/mK 

-2300  W/mK  [19] 

-790  W/mK  [20] 

(6,  6) 

3215  W/mK 

(8,  8) 

1813  W/mK 

250-350W/mK  [21] 

(10,  10) 

804.8  W/mK 

2 15-831  W/mK  [22] 

-2980  W/mK  [23] 

The  simulation  shows  that  the  thermal  conductivity  is  very  high  along  the  tube 
axis,  3325.8  W  m'  K'1  for  (5,5)  CNT  which  is  indicative  the  anisotropic  nature  of 
nanotubes.  This  value  is  comparable  to  those  obtained  by  Osman  et.  al,  whose  MD 
simulation  for  (5,  5)  SWNT  yielded  a  value  of  2300  W  nfi'K"1  [19].  The  trend  obtained 
from  this  simulation  shows  steady  decline  in  thennal  conductivity  as  the  diameter  of  the 
nanotube  increases.  A  similar  observance  is  made  by  Osman  and  Srivastava  for  (5,5)  and 
(10,  10)  nanotubes.  They  found  (5,5)  nanotubes  to  provide  the  highest  thennal 
conductivity  at  room  temperature,  as  compared  to  (10,  10)  and  (15,15)  nanotubes.  It  was 
suggested  that  the  probability  of  Umklapp  processes  are  higher  in  SWNTs  of  smaller 
diameters  as  compared  to  tubes  of  larger  diameters  [19].  The  Umklapp  processes 
describe  the  resistive  phonon-phonon  scattering  process  believed  to  limit  thennal 
conductivity  in  crystalline  materials. 

Unfortunately  there  is  little  agreement  on  the  specific  thermal  conductivity  of  a 
SWNT  as  indicated  by  table  2.  Typical  values  obtained  from  MD  simulations  have 
ranged  between  2000  and  6000  W/mK.  Further  investigation  is  required  to  explain  the 
trends  found  in  this  simulation. 
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B.  THERMAL  CONDUCTIVITY  OF  TWO  CARBON  NANOTUBES 

The  non-bonded  equilibrium  distance  between  equivalent  adjacent  SWNTs  was 
determined  from  the  Lennard-Jones  potential  function.  Figure  7  shows  the  LJ  potential 
curve  for  a  pair  of  armchair  (6,  6)  carbon  nanotubes.  The  equilibrium  distance  values 
obtained  are  shown  in  Table  3. 
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Figure  7.  Non-bonded  equilibrium  distance  between  two  (6,  6)  CNTs  oriented  at 

various  respective  angles. 


Table  3.  Lennard-Jones  Potential  Results  for  (6,  6)  SWNT 


Relative  Angle 

Equilibrium  Distance  (A) 

Minimum  Energy  (Kcal/mol) 

0° 

3.1 

7.228  eV 

15° 

2.86 

5.271  eV 

30° 

2.86 

3.052  eV 

45° 

2.84 

2.167  eV 

o 

O 

2.86 

1.792  eV 

75° 

2.86 

1.606  eV 

90° 

2.86 

1.538  eV 

26 


By  plotting  the  minimum  energy  distance  as  a  function  of  relative  angle  (Figure 
8),  it  is  apparent  that  beyond  15°,  the  relative  angle  does  not  have  a  significant  effect  on 
minimum  energy  position.  This  is  to  be  expected  because  a  cutoff  distance  of  r  =  2.5cr 
was  implemented  to  decrease  the  computation  time  of  the  LJ  potential.  As  the  angle  is 
increased  the  number  of  atoms  interacting  becomes  approximately  equal.  Although  the 
above  arguments  are  considered  for  only  two  nanotubes,  similar  arguments  apply  to  the 
entire  network  of  non-bonded  nano  tubes. 
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Figure  8.  Minimal  potential  energy  separation  distance  as  a  function  of  orientation 

angle,  6 


Once  the  minimum  equilibrium  distance  was  obtained,  the  thermal  conductivity 
of  paired  SWNTs  is  obtained  by  molecular  dynamics  simulation. 


Table  4.  Thermal  Conductivity  of  Two  (6,  6)  SWNTs  at  Various  Relative  Angles 


Angle 

Thermal  Conductivity 

0° 

15  W/mK 

15° 

20  W/mK 

75° 

15  W/mK 

VO 

o 

o 

10  W/mK 
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C.  EFFECTIVE  THERMAL  CONDUCTIVITY  OF  RANDOMLY  ORIENTED 
CARBON  NANOTUBES 

The  longitudinal  and  transverse  thennal  conductivity  values  obtained  from  the 
molecular  dynamic  simulations  conducted  previously  were  used  as  material  properties  in 
the  finite  element  analysis  of  the  effective  thermal  conductivity  of  randomly  dispersed 
SWNTs  in  a  unit  area. 


Table  5.  Material,  Geometric  Properties  and  Boundary  Conditions  for  FEA 


Conductivity  in  Axial  Direction 

3000  W/mK 

Conductivity  in  Radial  Direction 

20  W/mK 

Radius 

0.01 

Spacing 

0.01 

Tolerance 

0.01-0.1 

Boundary  Temp.  1 

100K 

Boundary  Temp  2 

200K 

In  Figure  9,  the  effective  thermal  conductivity  of  the  randomly  dispersed 
nanotubes  increases  nonlinearly  with  nanotube  loading.  Figures  10  and  11  shows  a 
steady  increase  in  the  effective  thennal  conductivity  in  both  the  horizontal  and  vertical 
directions  of  the  unit  square  as  the  length  and  total  volume  of  CNTs  increase.  In  addition 
Figure  9  shows  an  increase  in  effective  thennal  conductivity  with  increasing  tube  length. 
Finally,  Figure  12  shows  a  steady  increase  in  effective  thennal  conductivity  as  the 
minimal  number  of  CNTs  linking  the  boundary  is  increased.  Collectively  these  results 
show  that  the  heat  flux  between  boundaries  increases  with  tube  length  and  volume 
fraction.  As  the  volume  fraction  of  nanotubes  increases,  the  number  of  interactions 
increases.  In  turn  the  effective  conductivity  increases  due  to  more  efficient  percolation. 
The  aforementioned  behavior  implies  that  thennal  conductivity  enhancement  in  nanotube 
suspensions  will  largely  depend  on  length  and  concentration  of  nano  tubes.  The  findings 
presented  in  this  simulation  are  in  agreement  with  experimental  results  published  by 
Biercuk  et  al  which  show  nonlinear  increase  in  thermal  conductivity  with  an  increasing 
weight  percentage  of  SWNT  loaded  epoxy  [5]. 
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Figure  9.  Enhancement  in  effective  thermal  conductivity  as  a  function  of  the  number  of 

CNTs 
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Figure  10.  Enhancement  of  effective  thermal  conductivity  in  the  horizontal  direction  as  a 
function  of  the  total  volume  of  CNTs  with  different  CNT  lengths 
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Figure  11.  Enhancement  of  effective  thermal  conductivity  in  the  vertical  direction  as  a 
function  of  the  total  volume  of  CNTs  with  different  CNT  lengths 


Min.  No  of  CNTs  Connecting  Boundary 


Figure  12.  Enhancement  of  effective  thermal  conductivity  in  the  horizontal  and  vertical 
direction  as  a  function  of  the  no.  of  CNTs  linking  the  boundary. 
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V.  CONCLUSION/RECOMMENDATION 


The  heat  conduction  coefficient  of  a  single-walled  carbon  nanotube  (SWNT)  was 
simulated  by  the  molecular  dynamics  method  with  the  Tersoff-Brenner  bond  order 
potential.  The  results  obtained  from  the  simulation  were  then  used  as  material  properties 
in  the  finite  element  model  of  the  effective  thermal  conductivity  of  randomly  oriented 
carbon  nanotubes.  This  method  is  based  on  a  system  of  one-dimensional  percolating  rods 
that  connect  two  reservoirs  at  different  temperatures.  The  significant  findings  of  this 
work  were: 

•  A  multiscale  technique  is  useful  to  predict  the  effective  thermal 
conductivity  of  randomly  dispersed  CNTs. 

•  The  equilibrium  distance  between  two  adjacent  CNTs  remains  almost  the 
same  as  the  relative  angle  varies. 

•  For  the  same  size  of  CNTs,  the  effective  thermal  conductivity  increases 
nonlinearly  with  the  increasing  number  of  CNTs. 

•  The  effective  thennal  conductivity  is  significantly  greater  for  longer  CNTs 
as  compared  to  shorter  CNTs  with  the  same  total  volume  fraction. 

The  general  trends  resulting  from  these  simulations  indicate  significant  thermal 
conductivity  enhancement  when  the  tube  length,  volume  fraction  and  number  of  CNTs 
linking  the  boundary  is  increased.  These  results  support  the  very  important  implication 
of  thermal  enhancement  in  nanotube-based  composite  materials  as  well  as  nano-devices. 

A  few  recommendations  can  be  made  for  future  studies  based  on  the  methods 
used  and  results  obtained  in  this  research. 

•  In  reality,  carbon  nanotubes  have  defects  and  irregularities  that  decrease 
the  thennal  conductivity.  Therefore  vacancy  concentrations  should  be 
incorporated  in  the  simulation  to  determine  its  affect  on  thermal 
conductivity. 

•  Research  has  shown  that  magnetically  aligned  nanotubes  in  composite 
materials  further  enhance  the  thermal  conductivity.  The  finite  element 
model  used  here  can  be  adjusted  to  investigate  these  findings  further. 
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APPENDIX 


CNT  MODELS 


Figure  13.  Relative  Angle  =  0°  Figure  16.  Relative  Angle  =  45 
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Figure  15.  Relative  Angle  =  30° 


Figure  18.  Relative  Angle  =  75 
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